Power-laws f(R) theories are cosmologically unacceptable 
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In a recent paper [lj we have shown that f(R) = R + pR n modified gravity dark energy models 
are not cosmologically viable because during the matter era that precedes the accelerated stage the 
cosmic expansion is given by a ~ t 1 ^ 2 rather than a ~ t 2 ^ 3 , where a is a scale factor and t is the 
cosmic time. A recent work by Capozziello et al. 0] criticised our results presenting some apparent 
counter-examples to our claim in f(R) = pR n models. We show here that those particular R" models 
can produce an expansion as a ~ t 2 ^ 3 but this does not connect to a late-time acceleration. Hence, 
though acceptable f(R) dark energy models may exist, the R n models presented in Capozziello et 
al. are not cosmologically viable, confirming our previous results in Ref. [J]. 



Among the various interesting possibilities invoked in 
order to explain a late-time accelerated expansion, f(R) 
modified gravity dark energy (DE) models (R is the Ricci 
scalar) have attracted a lot of attention. 

However, we found recently in Ref. [l[ that for a large 
class of f(R) DE models, including R n models, the usual 
power-law stage a(t) oc t 2 l z preceding the late-time ac- 
celerated expansion is replaced by a power-law behaviour 
a(t) oc t 1 / 2 . Such an evolution is clearly inconsistent with 
observations, e.g. the distance to last scattering as mea- 
sured by CMB acoustic peaks. Hence a viable cosmic 
expansion history seems to be a powerful constraint on 
such models. The results in Ref. [l| have been discussed 
in greater detail and largely expanded in Ref. [3[ . These 
cosmological difficulties of f(R) models are complemen- 
tary to those that arise from local gravity constraints. 

However, the claim of [l| was recently criticised by 
Capozziello et al. (CNOT) 0. This short paper is de- 
voted to addressing explicitly this criticism. 

It is necessary to begin with some clarifications. First, 
it is clear that f(R) gravity models can be perfectly vi- 
able in different contexts. Maybe the best example is 
Starobinsky's model, f(R) = R + pR 2 [4|], which has 
been the first internally consistent inflationary model. 
This Lagrangian produces an accelerated stage preced- 
ing the usual radiation and matter stages. A late-time 
acceleration in this model requires a positive cosmologi- 
cal constant (or some other form of dark energy) in which 
case the late-time acceleration is not due to the R 2 term. 

Second, it is important to clarify an issue raised in 
CNOT concerning the validity of the conformal trans- 
formation we used in We checked all our results 
numerically (and where possible also analytically) both 
in Jordan frame (JF) and Einstein frame (EF), always 
considering the former as the physical frame (i.e. the 
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frame in which matter is conserved with an energy den- 
sity p m oc a~ 3 ). So the power-law behaviour a ~ i 1 / 2 
found in JF is not an artifact of the conformal trans- 
formation. It is in fact the same solution of the original 
Brans-Dicke paper @ in 1961 (Eq. 60) with uj = (equiv- 
alent to = 1/2 in the notation of [l() and corresponds 
to solutions found in JF also in other R n papers such as 
Ref. @. 

The main criticism in CNOT is that it is possible to 
have a stage with a(t) oc < 2 / 3 followed by a DE dominated 
phase for some f(R) models and even for the power-law 
case f(R) = pR n (notice that we changed the sign of 
n with respect to our paper [jj in order to match the 
choice in CNOT) and several examples are suggested. We 
address here the viability of the R n models suggested in 
CNOT. 

The f(R) gravity action in the JF is given by 



S = 



1 



2 f{R) 



(i) 



where k 2 = 8irG while G is the gravitational constant, 
and £ m and £ ra( j are the Lagrangian densities of dust-like 
matter and radiation respectively. In the flat Friedmann 
Robertson Walker metric with a scale factor a, we get 
the following equations 

3FH 2 = K 2 (p m + p lad ) + l-(FR- f)-3HF, (2) 



-2FH = k z [p 



gPrad 



F-HF. 



(3) 



where F = df/dR and H = d/a with a dot being the 
derivative in terms of cosmic time t. Note that the energy 
densities p m and p ra d satisfy the conservation equations 
p m + 3Hp m = and p rad + 4i7p rad = 0, respectively. 
We shall introduce the following quantities 
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Then for the power-law models f(R) — (iR 71 , we obtain 



diV 
dx 2 
dF 



-1 + x\ — 3x 2 + nx 2 (\ + xi) + x 3 , (5) 
X1X2 + x 2 (xi + 2nx 2 + 4) , (6) 



n — 1 



^ = (xi - 2nx 2 )x 3 , 
where TV = In a. Wc also define 



fin 



3FH 2 



1 — Xi 



(1 - n)x 2 - x 3 . 



(7) 



(8) 



For general f(R) dark energy models one needs an addi- 
tional equation to close the system [3[. Among the fixed 
points which exist for the above three dimensional sys- 
tem, the following three types of solutions are important 
for our discussion (we assume no radiation, i.e. x 3 = 
in the following, unless otherwise stated). 

(i) Solution A ("curvature-dominated solution"): 
This corresponds to the fixed point 



(£1,2:2) 



2(2 - n) 



5 - An 



In - 1 ' (l-2n)(l-n) / ' 
The evolution of the scale factor is given by 

(1 - 2n)(l - n) 



fim = . 



a oc t aA , a a 



2-n 



(9) 



(10) 



It is an exact solution in the absence of dust, and an 
asymptotic solution in the presence of dust. The latter 
was originally used to give rise to a late-time accelera- 
tion for negative n ("curvature dominated late-time at- 
tractor" ) [3 Hj] • When a a < the expanding solution is 
given by a oc (t s — t) aA , which corresponds to a phantom 
solution. 

(ii) Solution B ("scaling solution"): 
This corresponds to the fixed point at which the energy 
fraction of the matter does not vanish: 



(xi,x 2 ) 



3(n-l) 



(An - 3) 



n 2n 2 
8n 2 + 13n - 3 
27? " 



(11) 



The evolution of the scale factor is given by 

a oc f B , a B = 2n/3. (12) 

(iii) Solution C ("ci matter dominated epoch" [10]): 
This corresponds to the fixed point 



(x u x 2 ) = (-1,0) , fi m = 2. 



(13) 



This stage is the so-called ^-matter-dominated era 
(cAMDE) 10] with scale factor evolution 



for any n. 

It was shown in [l[ that for all n the </>MDE replaces 
the usual matter era prior to the late-time acceleration 
driven by the point A. CNOT instead pointed out that 
it is possible to use either solution A or B in order to 
have a standard matter era (a oc t 2 / 3 ) followed by an 
accelerated expansion. Clearly, two possibilities arise: 
either the universe goes from A to B or vice versa. In the 
first case the solution A has to behave as a matter era 
(a A = 2/3), and therefore n = -0.129 orn = 1.295. In 
the second case we require the condition as = 2/3, which 
corresponds to n = 1. Hence the three possible "counter 
examples" suggested by CNOT are: n — —0.129, n = 
1.295 and n = 1. Now we shall investigate whether these 
cases really provide a viable cosmological evolution. 

Let us first analyse the stability of the solutions A and 
B. Neglecting radiation and considering linear perturba- 
tions around the fixed points (x%, x 2 ) along the line pre- 
sented in Ref. [ll| , we obtain the following eigenvalues of 
the corresponding Jacobian matrix evaluated at the each 
fixed point: 



An 



Hi 



1 - n 

for the solution A, and 



M2 



8n 2 - 13n + 3 
"(l-2n)(l-n) 



(15) 



M± = 



aott ac , a c = l/2, 



(14) 



3(1 - n) ± - n)(-256n 3 + 608n 2 - A17n + 81) 

An(n — 1) 

(16) 

for the solution B. 

This shows immediately that the case n — —0.129 (and 
values in the vicinity) is excluded because the point A 
is then stable (u\ < 0, u 2 < 0) and, once reached, it 
will never give way to a late-time acceleration. In other 
words, the transition from A to B is impossible in this 
case. When n = 1.295, A is a saddle point (/xi < 0, \i 2 > 
0) and B is a stable spiral. Hence the the transition from 
A to B is possible, but with this value of n (and values 
in the vicinity) the point B is not accelerated, since then 
glj3 = 0.863. It is also important to note that the point 
A corresponds to a solution without matter (fi m = 0), 
so this would be a "matter era" without matter, which 
is clearly not acceptable as well. This leaves as the only 
possibility n — 1 and a transition from B to A. 

From Eq. ^ the point A disappears for n — 1, which 
means that the transition from B to A is not possible. 
As this case merely corresponds to Einstein gravity, it is 
obvious that one gets the required behaviour a oc t 2 ! 3 in 
the dust-dominated era. However there is no mechanism 
left for the generation of a late-time acceleration unless 
some additional DE component is introduced, which is 
what modified gravity DE models are supposed to avoid. 

So we conclude from the discussion above that the so- 
lutions suggested in CNOT do not lead to a a oc t 2 ^ 3 
behaviour followed by an accelerated expansion. 

Still it may be interesting to consider the scenario with 
n close to 1, for which a oc t p with p w 2/3, instead of 
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exactly 2/3. Let us study the case with n in the con- 
servative range 0.75 < n < 1.25, which corresponds to 
power-law exponents 1/2 < p < 5/6. Since the point B 
is a stable spiral for 1 < n < 1.327, transition from a 
decelerated matter era to an accelerated one is impos- 
sible also in this case. For 0.713 < n < 1 the point B 
is a saddle, so a transition is indeed possible. For these 
values the point A corresponds to a stable node with an 
effective phantom equation of state 



(n 7^ 0) models since the reconstruction attempted in 
CNOT is not for f(R) models of this type. 



w e f[ = -1 



2(2- 



3(1 - 2n)(l - n) 



< -7.6. 



(17) 



Meanwhile the third point C, the </>MDE, is a stable point 
as well with an effective equation of state w e s = 1/3. 
Then the trajectories leaving the point B are attracted 
either by A or C. However the final accelerated point A is 
a strongly phantom one as given in Eq. (JTTJ) . In addition 
the more one tries to obtain the exact matter era (n — ► 1) 
the more phantom the final acceleration is (w e g — > — oo 
as n — > — 1). Hence from this point of view these models 
are cosmologically unacceptable. 

To complete our proof, we have run our numerical code 
to investigate the evolution of the system in the space 
(x\,X2)- Without including radiation the final attractor 
is in fact either A or C depending upon initial condi- 
tions. However trajectories which are attracted to B first 
and then finally approach the point A are restricted in a 
narrow region of phase space (see Fig. [T|) . Moreover the 
duration of the (quasi) matter epoch gets shorter as we 
choose the values of n closer to 1 (which is in fact required 
to obtain the matter phase). This is associated with the 
fact that the eigenvalues of the matrix of perturbations 
diverge in the limit n — > 1 [|[ . When we start from realis- 
tic initial conditions <C 1, |ac2 j "C 1) with inclusion of 
radiation, the solutions directly approach the fixed point 
A or C without passing the vicinity of the point B. In 
other words we have not found any trajectories in which 
the radiation era is followed by matter and final accel- 
erated eras. Therefore, our numerical analysis excludes 
also the range 0.713 < n < lasa viable cosmological 
model. 

In CNOT the possibility is also mentioned of recon- 
structing the theory from observations (in particular from 
the function H(z) where z is the redshift), in analogy to 
the reconstruction of scalar-tensor DE models [H, EH]. 
Clearly we expect that many f(R) DE models can be 
successfully reconstructed at low redshifts for any H(z) 
corresponding to late-time accelerated expansion. How- 
ever, nothing guarantees that an H{z) corresponding to 
a conventional cosmic expansion at high z leads to an ac- 
ceptable form of f(R). Sometimes a singular behaviour 
arises already at very low redshifts [HI, E3 • The proce- 
dure of reconstruction does not guarantee either the sta- 
bility of the solution on higher redshifts. Finally, the par- 
ticular reconstructed model proposed in CNOT does not 
contradict our claim on f(R) = fiR n or f(R) — R + \xR n 



f(R)=R 
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Figure 1: Phase space plot in the plane (srijXa) in Poincare 
coordinates for the model f(R) = /iR ' 9 in the absence of ra- 
diation. The dotted lines represent trajectories at the early 
stage, whereas the continuous lines correspond to those at the 
final stage. We show the fixed points A, B and C as circles. 
The final attractor is either the curvature-dominated solution 
A or the 0MDE point C. In the absence of radiation it is pos- 
sible to have the sequence from the (quasi) matter point B to 
the phantom point A with restricted initial conditions. As n 
gets closer to 1, however, the matter era becomes very short 
and the effective equation of state for the point A diverges 
(w e s — > — oo). If we take into account radiation (correspond- 
ing to initial conditions \xi\ <C 1, \x2\ -C 1 and xs ~ 1) the 
solutions do not approach the (quasi) matter point B before 
reaching the final attractor, which means the matter epoch is 
absent. 



The question of the cosmological viability of f(R) mod- 
els is a very interesting one. In Ref. [3j we have spelled 
out the conditions on the forms of f(R) in order to sat- 
isfy the basic cosmological requirements of a standard 
matter era and a late-time accelerated attractor. The 
specific f(R) = fiR n examples suggested in CNOT do 
not pass these criteria and therefore confirm, rather than 
contradict, our claim. 
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